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Abstract. In pure anisotropic or layered superconductors thermal fluctuations induce a van der Waals 
attraction between flux lines. This attraction together with the entropic repulsion has interesting conse- 
quences for the low field phase diagram; in particular, a first order transition from the Meissner phase to 
the mixed state is induced. We introduce a new variational approach that allows for the calculation of the 
effective free energy of the flux line lattice on the scale of the mean flux line distance a, which is based on 
an expansion of the free energy around the regular triangular Abrikosov lattice. Using this technique, the 
low field phase diagram of these materials may be explored. The results of this technique are compared 
with a recent functional RG treatment of the same system. 

PACS. 74.60.Ec Mixed state, critical fields, and surface sheath - 74.72.IIs Bi-based cuprates 



! 1 Introduction 

■ The physics of flux lines in high-Tc superconductors at 
' low magnetic fields is dominated by the competition be- 
tween the bare repulsion between flux lines and an en- 
tropic repulsion due to thermal fluctuations. This results 
, in a first order melting transition from the Abrikosov lat- 
' tice [|] to a liquid phase close to the lower critical field 
I iJcn as has been predicted some time ago by Nelson 
• For mean flux line distances, a, larger than the London 
' penetration depth A, the bare interaction decays expo- 
, nentially with a, while the entropic repulsion decays al- 
' gebraically ~ {X^ / Lto)'^ . Here, the thermal length scale 
Lt = <Pl/{16Tr^T) « 2cmK/r denotes the length of 
1 an isolated flux line segment that shows a thermal mean 
' square displacement of the order of A [||, and = hc/2e 
. is the flux quantum carried by each flux line. Close to the 
' transition, the entropic repulsion dominates over the bare 
I interaction, leading to a magnetic induction B that van- 
ishes linearly with the reduced field strength h = {H — 

In strongly anisotropic or layered superconductors an 
additional interaction between fiux lines has been found 
recently by Blatter and Geshkenbein (BG) Short scale 
fluctuations on the scale of pancake vortices lead to an 
attractive van der Waals (vdW) interaction. (This is in 
fact also the origin of the long range attraction of a flux 
line to the surface ||^.) For flux lines separated by a dis- 
tance R the strength of this interaction is of the order 
-X^/idLrR'^) ioi X < R < d/e and ~ -XyieLrR^) 
for d/e < R < X/e. Here, = ra/M <C 1 denotes 



the anisotropy of the material with m and M the ef- 
fective masses parallel and perpendicular to the Cu02 
plane, respectively, and d the interlayer spacing. Apart 
from this thermally induced attraction, also frozen-in dis- 
order in impure superconductors induces an attraction be- 
tween flux lines with the same dependence on the dis- 
tance R (Hj^J . This disorder induced attraction dominates 
at very low temperature T <C Tdis, where Tdis depends on 
the disorder strength |3|. Here, we will focus on the oppo- 
site case T ^ Tdis where it is sufficent to consider thermal 
fluctuations. 

The competition among the bare, the entropic and the 
vdW interactions leads to an interesting phase diagram 
at low B values. In particular, the vdW attraction can 
lead to an instability of the Abrikosov lattice in the dilute 
limit a ^ A, resulting in a first order transition between 
the Meissner and the mixed phase Q. 

In order to calculate the low B phase diagram of pure 
layered high-Tc superconductors, one has to calculate the 
Gibbs free energy density of this system on the scale of 
the mean flux line distance a, 

g{a-H,T) = f{a,T)-^-^K (1) 

which has to be minimized with respect to a with the 
external magnetic field H fixed. Here, f{a,T) is the free 
energy density, So — (^o/4vrA)^ — LtT/X^ is the basic 
energy scale, and n is the Ginzburg-Landau parameter. 
Note that the bare lower critical field (in the absence of 
thermal fiuctuations) is given by iJ"^ = A-kEo Iuk/^Qi and 
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that the magnetic flux is related to the flux hne distance 
via B ~ <?o/a^- 

It is the aim of the present paper to propose a new vari- 
ational approach to obtain this effective Gibbs free energy 
density, based on a perturbative expansion in small dis- 
placements of the flux lines around the regular triangular 
lattice. This expansion is justified only if we assume that 
the lattice is stable, i.e., that we are above the lower melt- 
ing line below which the flux lines form a liquid. The vari- 
ational technique allows for a check of this assumption by 
virtue of a Lindemann criterion where the self-consistently 
determined flux line fluctuations are compared with the 
mean flux line distance a. 

In order to motivate our efforts, let us shortly review 
the approaches that have been followed so far. In their 
first paper [Q, BG simply added the bare vdW energy, 
evaluated on the scale a, to the free energy which contains 
the bare and the entropic repulsion mentioned above. In 
this way, only contributions from the vdW energy on the 
length scale a are taken into account, which leads to a 
gross underestimation of its influence because of its rapid 
decay for i? > A. In |6|, in contrast, only the much more 
important contribution on the scale i?„iin was taken into 
account within a scaling approach, where i?min is the po- 
sition of the global minimum of the bare potential V{R). 
This scaling approach is elegant and simple and yields the 
correct qualitative picture, but does not allow for quan- 
titative predictions. More recently a functional renormal- 
ization group (RG) calculation has been adapted to the 
same problem |^]; in particular, the low field phase di- 
agram of Bi2Sr2CaCu20a; (BiSCCO) was calculated. In 
that approach, the bare flux line interaction is renormal- 
ized by thermal fluctuations on all scales between A and 
a, providing an effective free energy on the scale of a. The 
problem in this latter method is that the largest scale up 
to which fluctuations are taken into account is arbitrary to 
some extent; we will see below that this flaw is overcome 
in a natural way by the technique that we will develop 
below. 

The paper is organized as follows: In Section ^ the 
variational method is set up for a arbitrary bare flux line 
interactions. The method is applied to a flux line lattice 
in a layered high-Tc superconductor in Section ^ There, 
the flux line interaction is modeled by the superposition 
of the bare, short range repulsion and the long range van 
der Waals attraction, using physical parameters typical 
for BiSCCO. The data for the free energy obtained there 
are used for the minimization of the Gibbs free energy 
density in Section ^. We shortly review the functional RG 
approach in Section ^ and compare the two methods. Con- 
clusions are drawn in Section pi 



the mean distance between the FLs. Since in equilibrium 
the FLs form a triangular lattice we choose 

a„ = a(cos(0„)e^ -f sin((/)„)ej,), '^"^Y2^"'3' ^"^^ 

For n — 1, . . . , 6 these a„ are the distance vectors R — R' 
for all nearest neighbors R' of R on a triangular lattice. 
With this definition of the vectors a„ the equivalence of 
the X and the y direction is ensured. This results in the 
simple form for the matrix K°'^{ci) defined below. 

We write the actual position of line (n, m) at height z 

as 

The interaction between the vortices {n,m) and {n' ,m') 
is described by the function F(r(z)) that depends on the 
distance at equal height z. We have made the usual as- 
sumption that the vortex coordinates vary slowly with z, 
so that the FL interactions are well approximated by a 
potential which is local in z] this approximation is only 
strictly valid when the vortices are parallel to z. 
The system is governed by the Hamiltonian 

"'0 n=lm-l 

+ '^^ ^ V{Vn^rn{z) -Vn',m'{z))- (4) 

^ (n,m)7^(n', 771') 

This description of a FL lattice applies in the limit of fluc- 
tuations on scales larger than A, with a non-dispersive line 
stiffness k = £i{kz <C 1/A); on smaller scales, the FL stiff- 
ness £i(kz) in anisotropic superconductors is in fact highly 
dispersive [||. The FL interaction V{R) will be specified in 
section |3[ the following treatment makes no assumptions 
about its specific form. 

Applying periodic boundary conditions in all 3 spatial 



and u. 



directions (in particular u„_|_jv,7« = t 
U7i.,7i), Eq. (0) can be rewritten in Fourier space as 



n,777+Af 







dz J2 d\ 



X 



(5) 



q {n.m)^{n' ,m') 

The wave vector qz takes the discrete values 0, 27r/L|| , . . . , yl. 
The ultra violet cutoff A has to be introduced explicitly 
in order to avoid a divergence of the free energy. For the 
discrete set of — {qx,qy) vectors that are summed over, 
we choose 



2 Variation of the Free Energy 

Consider a system with N times M flux lines (FLs) in 
D = 3 dimensions, confined to a box of base area and 
height L\\ . The equilibrium position of line number (n, m) 
is given by R7i,rn = nsii + ma.2, where a = |ai| = |a2| is 



2 2n 
V3 a 



/ — cos(7r/12) sin(7r/12) 



(6) 



M 



cos(7r/12) 



N 



sin(7r/12) 



with n = 0, . . . , TV— 1 and m = 0, . . . , Af — 1. This is equiv- 
alent to choosing the vectors from the first Brillouin 
zone (in the 2 dimensional plane perpendicular to the field 
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Fig. 1. First Brillouin zone and summation area in q^-space. 
The vectors q with (n,m) = (iV, 0), (0, M) and {N,M) as 
defined in Eq. (pi) are reciprocal lattice vectors. 



direction). The summation area in Fourier space, together 
with the first Brillouin zone, is iUustrated in Fig. |l|. 

The Fourier transform of the fluctuating field u„^„i(z) 
is defined as 



and the Fourier transform of the potential V{r) as Vk 
/ (fxV{x)e-'^''. 

Now, we define a Gaussian variational Hamiltonian 

= ?E^"''(q)""(q)^'5(-q)' (8) 



with summation over the Cartesian components a, P — 
1,2, and rewrite the original Hamiltonian as 

H = Ho + iH- Ho). (9) 

Let (. ..)o denote thermal averages with respect to Hq. 
Clearly, 



where K{q) is the Matrix inverse of K{q), i.e.. 
The correlation function K{q) is a 2x2 matrix: 



(10) 



K(q) 



K*o{q) Koiq) 



(11) 



While the diagonal terms X£)(q) are real, the off-diagonal 
part Ko{q) may in principle take on complex values. In 
a straightforward but tedious calculation it can be shown 
however that Ko{q) is actually real, too. For complete- 
ness, we give the explicit expression for the matrix inverse 
^(q): 

■ft^D(o)(q) 



K 



D{0) 



(q) = +(-) 



Kl{q)-Kl{q) 



(12) 



As stated in the introduction, our aim is to minimize 
the Gibbs free energy density g{a] H, T) as defined in Eq. (||) 
as a function of the mean FL separation a. Let us first 
consider the free energy F = F{a,T). We define an ap- 
proximate free energy F* that is obtained by a first order 
perturbation expansion around Fq, the free energy associ- 
ated with Hq, 



F* =Fn 



(H-H 



0/0- 



(13) 



It is well known Q that for any choice of Hq, F* is an 
upper bound on the true free energy F. Therefore min- 
imization of F* with respect to the parameters defining 
Ho, Koiq) and Ko{q), yields the best approximation of 
its kind, namely Hq corresponding to the optimal param- 
eters is the best fluctuating lattice Hamiltonian that may 
be used to describe our physical system. The free energy 
Fo is given by 

Fo = -r'ln j PHe-'^^" ^ i_^lndet(/3if(q)). 

(14) 

Next, the expectation value {H — Ho)o shall be calcu- 
lated. To proceed, it is useful to rewrite the last term in 
the interaction part of (^ using the Fourier representation 
(0). Then, the averaging can be immediately carried out, 
leading to 



(7) {H - Ho)o = f E Tr [X(q) ] - L.NM^ksT 



L» NM 



e'^ X 



(15) 



We restrict the sum to a summation over nearest 
neighbors, which is a good approximation in the dilute 
limit we are interested in; in principle, one can also include 
next-nearest neighbors, etc. For the triangular lattice, it is 
convenient to choose the 6 vectors a„ from Eq. (H), with 
n = 1, . . . , 6. 

In the following, we choose the temperature ksT — 
(3^^ as the basic energy scale, and introduce the parame- 
ters K, = the reduced potential u(r) = I3V{y) and the 
reduced Hamiltonian Ti — f3H. Finally, we replace f3K{q) 
by i^(q). 

The variational parameters Kd (q) and Kq (q) are de- 
termined by minimizing F*: 



dF* 



dKD[o){q) 



0. 



(16) 



This leads to 



Koiq) - kql - 2^sin2(q^ • a/2) A, (17) 

a 

i^o(q) =-2^sin2(q^ .a/2)i?a (18) 
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where we have defined 



1 
1 



Aa = ^ / k2w(k)e-*'^e-^'=°^^»'''=" and (19) 



-i^ia-hkcM^'^kf, 



(20) 



By comparison with Eq. (p^, the 2x2 matrix Ma is defined 

as 



MNL, 



-^sin2(q^ .a/2)(w„(q)?/^(-q))o, (21) 



where we have plugged in the representation ( p^ of the 
matrix inverse K{q). Again, Ma,_D and Ma,o are the di- 
agonal and off-diagonal terms, respectively. Afa denotes 
the matrix inverse of Ma. Eqs. (|6|) to (^l]) serve as a 
self-consistent set of equations for K{cO- 



3 Van der Waals Interaction 

To be specific, we take the bare potential v{R) to be given 
by the superposition of the short range attractive and 
the long range repulsive interaction describing the direct 
vortex- vortex interaction in the extremely decoupled limit 
e ^ 0, 



v{R) = wo Ko{R/X) - a,d^(b{R/X) 



i?4 



(22) 



where vq = 2£o /ksT measures the amplitude of the direct 
interaction between flux lines, and ttvdw determines the 
strength of the thermal vdW attraction. (j){x) is a function 
that smoothly cuts off the power law part for i? < A, which 
we have defined as 




(' 



X < Xl 

Xi < X < X2 
X > X2 



(23) 

with Xl — 1 and X2 ~ 5. The choice of the cutoff function, 
as well as the actual values of Xi and X2, is to some extent 
arbitrary; a:i = 1 is however an obvious choice, and x^ 
has to be chosen such that the cutoff is not too sharp 
and, on the other hand, does not influence the form of the 
potential in the vicinity of the minimum for those values 
of flvdw that are physically meaningful . 

In order to make quantitative predictions for the low- 
flcld phase diagram of layered high-Tc superconductors, 
we identify the parameters introduced above with physi- 
cal parameters characterizing those systems: The elastic 
constant is k = eo/2kBT in the long wavelength regime 
Xkz <C 1 [D, and the amplitude of the vdW attraction 
(relative to wq) is given by avdw ~ kBT/{2eod\v?{'KX/d)), 
where d is the layer spacing 

Let us quantify these parameters for a speciflc highly 
anisotropic material, BiSCCO. This superconductor is char- 
acterized by the London penetration depth A ~ 2000 A, 



a Ginzburg-Landau parameter k « 100, a layer spacing 
d w 15 A and an anisotropy parameter e w 1/300 [||. To- 
gether with the thermal length Lt ~ 2 cm K/T, we flnd for 
our model parameters the values Xk w 10'°^ K/T, vq — 4k, 
and flvdw « 2 X 10"^ T/K. At T = 100 K, which is of 
the order of the critical temperature Tc, we finally have 
Xk w 10'^ and Cvdw ~ 2 x 10"'^. We want to feed the poten- 
tial (^2|) with K and vq discussed above into the variational 
procedure. The London penetration depth A is chosen to 
set the length scale in the direction perpendicular to the 
FLs, i.e., all lengths are measured in terms of A. 

Although the vdW amplitude avdw is determined by 
material parameters and the temperature as noted above, 
we will tune this parameter here in order to find whether 
a 'critical' value a*^^ exists where the phase transition 
from the Meissner phase to the mixed phase changes its 

to a first 



< a: 



vdw^ 



character from a second order (for a 
order transition (for a > a*^^). 

With these parameter values and the potentia l (22| ), 



we have solved the set of self-consistent equations (|l£ ) to 
( pl| ) numerically for different values of the mean spacing 
a with fixed FL number N x M and fixed vertical length 

In order to speed up computation, the sums over 
were replaced by integrals. 

First, we present the fluctuation matrices (^ as a 
function of the mean FL distance a. There are 6 fluctua- 
tion matrices corresponding to the 6 nearest neighbor vec- 
tors given by Eq. (||) . The special choice of these vectors 
a„, however, yields only two different matrices, namely 

M(^) = Afa„ for n = 1, 3, 4, 6 and 

M(2) ^ Afa„ for n = 2, 5. (24) 

With M*^^-' and M*^^-' each having two independent entries 
(the diagonal and off-diagonal terms Md and Mq, respec- 
tively), we have to solve for a total of 4 variables. The 
implementation made use of a multi-dimensional Newton- 
Raphson method |]l0j. Numerical solutions for these four 
observables are shown as a function of the mean spacing a 
in Fig. 1^. The data shown in this figure stems from calcu- 
lations with the vdW amplitude Ovdw — 10~^; the corre- 
sponding data for higher values of Ovdw is not significantly 
different. 

For values of a much smaller than the minimum posi- 
tion of the bare potential (which is at R^in ~ 27A for the 
vdW amplitude chosen here) fluctuations are strongly sup- 
pressed by the direct repulsion from the nearest neighbors. 
The fluctuations grow rapidly around a « 25A, because in 
that regime the strong repulsion is neutralized by the vdW 
attraction. For a large intermediate range, the fluctuations 
scale almost ^ a^. Note that in that plot, this quadratic 
scaling would correspond to a horizontal line. Actually, 
the slope is slightly negative, because the attraction con- 
tinuously decreases for larger distances. This tendency is 
most obvious in the plot of the fluctuations 



(u(0)2) = 



MNL« 



q 



^D(q) 



(25) 



as a function of a in Fig. For a > 150A, the size of the 
fluctuations saturates. This is due to the finite system size 
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Fig. 2. Fluctuations M^*^ divided by , as a function of the 



spacing a. The two curves M^j^^ and M)^' are hardly distin- 
guishable. The vdW amplitude has been set to Ovdw = 10~®, 
and the other model parameters to k = 10^, vq = A x 10^, 
N = M = 100, L|| = 5 X 10'', A/-K = 400. In the upper right 
corner of the figure, the function L\\ /2ir'^ka? has been added 
{dashed line) which gives the upper bound for as derived 

in Eq. @. 
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Fig. 3. The fluctuations ({u^(0)))^''^/a, as a function of a, 
from the same numerical data as in Fig. bl 



L\\ in the z-direction; in this Umit, the contribution of the 
potential energy to the free energy density is neghgible, 
leaving only the kinetic term. This leads to 

2tt kql ^ 27r2 K ' ^ 



277/L|| 



where we have assumed that the correlations between neigh- 
bors are small and that L\\A/2'k 3> 1. Finite size effects 
will hence occur when the mean vortex distance a becomes 
of the order of L\\ j^-K^k. In this regime, the off-diagonal 

terms Mq ^ and Mq^-* vanish due to the effective rota- 
tional invariance for each FL in this limit, where the FLs 
no longer interact with each other. 



Fig. 4. Gibbs free energy density g{a\ fi, T) for several val- 
ues of fj., with Ovdw = 10""". The data reflects a second order 
transition from the Meissner phase to the mixed phase. 



4 Minimization of the Gibbs Free 
Energy 

We will now use the effective free energy density /(a, T) = 
F{a,T)/L\\NMa?, where the free energy F{a,T) is taken 
to be the minimal value of F* , to calculate the Gibbs free 
energy density 



g{a;fi,T) = f{a,T) + fj./a'^ 



(27) 



which has to be minimized with respect to the mean FL 
distance a. For convenience, we have introduced the 'chem- 
ical potential' fi = —EohliiK, where h = {H — Hc-^) / Hc-^ is 
the deviation of the applied magnetic field from the lower 
critical field. 

In Figs. ^ and ||, the Gibbs free energy density g{a; fi, T) 
as defined in ( p7[ ) is plotted for two different values of Cvdw 
For these two values of the vdW amplitude, the transition 
from the Meissner state to the mixed phase has different 
characteristics: For Ovdw = 10~^, the transition is con- 
tinuous because the minimum position of the potential is 
continuously shifted to larger values with growing /i, cor- 
responding to a second order phase transition. 

For Cvdw — 10^**, on the other hand, the transition is 
discontinuous and hence first order, because for in the 
vicinity of /ic, two minima emerge: one at a finite value of 
a and the other at a = oo. The magnetization B ~ l/a^;^ 
as a function of H thus shows a first order transition from 
the Meissner phase to the mixed phase at a finite mag- 
netization By. Hence, the principal result from the study 
by BG Q and from the functional RG treatment Q is 
reproduced by the variational approach. 

The transition is illustrated in Fig. |[ where l/a^^^J^ii) 
is plotted as a function of /x—/ic. For each value of Ovdw, /^c 
has been determined individually such that g{a] fi > fic, T) 
has no minimum at finite a, while it does have such a min- 



imum for fi < He- For Cvdw < 



vdw 



the magnetic induc- 



tion is proportional to fi~ fic- This scaling arises from the 
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Fig. 5. Gibbs free energy density g{a;fj,,T), now with 
oivdw = 10~*. The figure reflects a first order transition from 
the Meissner phase to the mixed phase. 



competition between the bare repulsive interaction and 
the entropic repulsion alone. 

For the parameters chosen above, the 'critical' value 
of flvdw separating the two regimes is found to be a^^^ ~ 
2 X 10"^. As stated above, the physical value of the ampli- 
tude of the thermally induced van der Waals attraction is 



of the order 



10 , hence deep in the regime where 



the transition ~ according to the variational procedure - 
is first order. From the data in Fig. we read off that 
the mean separation between FLs just above the transi- 
tion between the Meissner phase and the mixed phase is 
Omin ~ 20A. For BiSCCO, this corresponds to a magnetic 
induction « 500 G/(ami„/A)2 w 1.2 G §. 

Now, let us face the question whether the basic as- 
sumption of the variational approach, i.e., the existence 
of a regular triangular lattice, is valid for these parame- 
ters. This can be checked using the self-consistently de- 
termined fluctuations by applying a Lindemann criterion 
||ll|,|2j. In particular, the elastic structure melts when the 
displacement between two neighbors becomes 



((U"(z) 



(28) 



where n and m denote two neighboring vortices, a is the 
lattice constant, and cl is the heuristically determined 
Lindemann number. In a recent numerical study on FL 
lattice melting |l^, this number has been found to be 
Cl ~ 0.25, in agreement with usual estimates for vortex 
lattice melting |8|. With the values represented in Fig. ||, 
we are hence above the lower melting line where the lattice 
is still stable ||l||. 

In principle, the variational technique can now be used 
to explore the whole low field phase diagram of layered 
superconductors, tuning the temperature T. We will re- 
strict ourselves here however to the data shown above, 
since we do not expect results for the phase diagram that 
differ qualitatively from those obtained by the functional 
renormalization group treatment to be discussed below. 



Fig. 6. 1/ a^in versus fi~ fic for several values of avdw In the 
lower left corner, a (dashed) line has been added that shows 
a linear scaling of the magnetic induction B as a function of 
fj,—fj,c, revealing the linear scaling of the two lower curves which 
correspond to the two lowest values of the vdW amplitude. 
The other curves approach a finite value for fi ^ fic, 

thus determining the van der Waals magnetic induction _B„ = 



Instead, we will compare the two techniques and discuss 
their relative strengths and flaws. 



5 Comparison with the functional RG 

As noted in the introduction, the problem considered here 
has recently been addressed by a functional renormaliza- 
tion group calculation J7|. In that approach, the FL inter- 
action potential, Eq. is renormalized by short scale 
thermal fluctuations, analogous to the first RG treatment 
of the purely repulsive short range FL interaction by Nel- 
son and Seung [ 12| . The starting point of this technique is 
the bare Hamiltonian with two interacting FLs. Due to 
the presence of neighbor vortices in a typical distance a, 
transverse fluctuations are assumed to be confined to this 
length scale. Hence, thermal fluctuations are integrated 
out by means of a momentum shell RG up to this scale, 
a, applying (nonlinear) recursion relations that are closely 
related to those that were established in the context of the 
wetting transition (Note that this argument applies 
to the system in a liquid state; in the solid state, the fluc- 
tuations are confined to a shorter scale < Clo, as we have 
seen above.) This procedure renormalizes the interaction 
potential, leading to an effective potential V'^^{R) on the 
scale a. 

For large a ;» A, the RG takes us into a region of weak 
coupling and high vortex densities |12 , where the effec- 
tive Gibbs free energy density simply reads g{a;fj,,T) = 
(|y^(a) + iJ,)/a^ . z = 6 is the number of nearest neigh- 
bors. Again, g(a;^,T) has to be minimized with respect 
to a in order to obtain the Gibbs free energy density as 
a function of fi and T, hence on this level the analysis of 
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Fig. 7. The effective interaction energy V^^ (a) which re- 
sults from the functional RG procedure when stopping the in- 
tegration on the scale a/rj. The data result from numerical 
functional RG calculations with a rescaling factor b = 1.2 and 
vo = 2 X 10^, corresponding to the temperature T = 100 K. 
The vdW amplitude is set to the values Ovdw = 10~*. 



the phase diagram is identical to that employed within the 
variational technique. This enables us to directly compare 
the results of the two approaches. 

Returning to the result for By obtained above for BiS- 
CCO at T w 100 K, one finds that this value is consid- 
erably higher than B^'^ w 0.2 G as calculated from the 
functional RG data . This discrepancy is consistent with 
the fact that the value for a^^^ determined above is much 



smaller than the value a*^^ « 2 x 10~^ which results from 
the RG calculations. 

The reason for the discrepancy between the two ap- 
proaches is given by the fact that they include fluctua- 
tions up to different maximal scales, as mentioned above. 
In particular, we have found here that fluctuations are 
confined to (u^) w (0.10 . . . 0.15)a. Hence, the FLs have 
a much smaller contact probability than assumed in the 
RG procedure, and the strong short range repulsion does 
consequently contribute less to the effective interaction 
energy, rendering the effective potential more attractive. 

This statement can be made more quantitative by ad- 
apting the data from the functional RG to the situation 
where the FLs form a lattice by confining the integration 
to scales < ^ a with a constant rj > 1. An example for the 

effective interaction V^-^ obtained in this way is shown in 
Fig. |. 

Indeed, the figure reveals that Vg^(a) exhibits a global 
minimum at a finite length scale a for large enough ry, while 
it lacks such a minimum for smaller 77 « 1. With growing 
77, the position of the minimum decreases, corresponding 
to larger values By. 

Let us compare the results which are shown in Figs. ^ 
and ^ from the variational and the RG approach, respec- 
tively, corresponding to the same temperature T — 100 K 
and the same vdW amplitude aydw = 10^"*. For this com- 
parison, we determine the appropriate coefficient i] from 



the magnitude (u(0)^) of the fluctuations, as measured in 
the variational calculation, with a set to the length amin 
that corresponds to the position of the free energy mini- 
mum. For the present set of parameters, amin ~ 25A (see 
Fig. The corresponding fluctuations have been deter- 
mined to be (u(0)2)i/2 ~ (0.1 ±0.02) a (see Fig. |). A rel- 
atively large uncertainty is attached to this value because 
(u(G)^) varies rapidly in this regime. Now, we would ex- 
pect the two techniques to yield comparable results when 
choosing 77 « 10 ± 2. Indeed, the data compare favorably 
well: V^^ (a) as determined from the functional RG with 
7/ = 8 exhibits a minimum at a/ A = 26 ± 2 (see Fig. 0), 
which is consistent with the value determined by the vari- 
ational procedure. Modified in this way, both techniques 
hence predict the same magnetic flux By at the transition 
from the Meissner phase to the mixed phase. 



6 Conclusion 

In layered high- Tc superconductors close to the lower melt- 
ing line _ffci where the mean flux line distance becomes 
larger than the London penetration depth A, a long ranged, 
attractive interaction between the vortices is induced both 
by thermal or by quenched-in disorder fluctuations which 
has interesting consequences for the phase diagram of these 
materials in the regime of low magnetic fields B. In par- 
ticular, this van der Waals attraction may lead to a first 
order transition from the Meissner phase to the mixed 
phase, both for pure and for disordered superconductors. 

Restricting ourselves to the pure case, we have estab- 
lished a variational technique that allows for the compu- 
tation of the van der Waals magnetic induction By , which 
is based on a self-consistent expansion around the regu- 
lar Abrikosov lattice. We have applied this method to the 
calculation of By in the case of BiSCCO, which is a typi- 
cal example of a strongly layered high- Tc superconductor. 
The results for the low field phase diagram of pure lay- 
ered superconductors are qualitatively and, to some ex- 
tent also quantitatively, consistent with the results from a 
functional RG approach. The latter is the more powerful 
tool in that it correctly reproduces logarithmic corrections 
to the leading scaling behavior. It is hence adequate for de- 
termining the true asymptotic behavior at the phase tran- 
sition. The variational technique, on the other hand, has 
the advantage that it provides an inherent control mech- 
anism as to where to stop the coarse graining process or 
- in other words - which is the scale on which the ef- 
fective interaction has to be considered. Furthermore, by 
self-consistently determining the flux line fluctuations, the 
results from the variational method can be used to check 
a posteriori, using a Lindemann criterion, whether an ex- 
pansion around the regular Abrikosov lattice is justified 
for a given set of physical parameters - which is the basic 
assumption of the variational calculation. Indeed, this has 
to be checked carefully since we are working in a part of 
the phase diagram where the lower lattice melting line Bm 
and the van der Waals magnetization By are of the same 
order. 
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The variational approach may straightforwardly be ex- 
tended to the situation with quenched-in impurities us- 
ing the replica trick. An advance in this direction could 
be guided by variational studies of flux lines lattices in 
high- Tc superconductors in the presence of weak disorder 

We thank Jan Kierfeld and Thorsten Emig for helpful discus- 
sions and the German Israeli Foundation (GIF) for financial 
support. 
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